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We study the spin-1/2 two and three dimensional Orbital Compass Models relevant to the problem of orbital 
ordering in transition metal oxides. We show that these systems display self-dualities and novel (gauge-like) 
discrete sliding symmetries. An important and surprising consequence is that these models are dual to (seem- 
ingly unrelated) recently studied models of p + ip superconducting arrays. The duality transformations are 
constructed by means of a path-integral representation in discretized imaginary time and considering its spa- 
tial reflection symmetries and space-time discrete rotations, we obtain, in a transparent unified geometrical way, 
several dualities. We also introduce an alternative construction of the duality transformations using operator 
identities. We discuss the consequences of these dualities for the order parameters and phase transitions of the 
orbital compass model and its generalizations, and apply these ideas to a number of related systems. 

PACS numbers: 71.27.+a, 71.28.+d, 77.80. -e 



I. INTRODUCTION 

Orbital compass models offer a simple and qualitative de- 
scription of the ordering of orbital degrees of freedom in a 
number of complex oxides such as the titanates (J. The de- 
grees of freedom of these models describe the spatial orien- 
tation of the orbital degrees of freedom. Jahn-Teller effects 
lead to anisotropic orbital compass like interactions amongst 
the orbitals. When combined with the spin degrees of free- 
dom, to which the orbitals are coupled via super-exchange in 
these systems as well as by spin-orbit interactions, they 
lead to complex phase diagrams with phases that involve both 
spin and orbital ordering (and disorder) to various degrees. 
Indeed, these systems offer an interesting laboratory for the 
study of interesting anisotropic quantum nematic phases, with 
and without spin order, and are a simple example of electronic 
liquid crystal phases 01 0- 

Orbital compass models also exhibit unusual and so far not 
well studied symmetries which play a big role in their phys- 
ical properties. In the current article, we elucidate the dis- 
crete "sliding" gauge-like symmetries present in the two and 
three dimensional orbital compass models. In two dimen- 
sions, these symmetries involve flipping the orbital degrees 
of freedom simultaneously along a single row or column of 
the lattice. These discrete symmetry transformations stand in- 
between the global symmetries familiar from spin systems and 
the local symmetries of gauge theories. Although these are 
not truly gauge symmetries in the sense that they affect the 
boundary conditions, they are softer than the familiar global 
symmetries. In fact, for reasons discussed elsewhere 01 these 
discrete sliding symmetries are alike gauge symmetries in that 
they cannot be spontaneously broken. A direct consequence 
of the existence of these discrete sliding symmetries is that 
their natural order parameters are nematic, which are invari- 
ant under discrete sliding symmetries. Here we give an ex- 
plicit construction of the nematic order parameters and poten- 
tial physical consequences are discussed. 

Gauge-like symmetries appear in a number of condensed 



matter systems. Exact local gauge symmetries are pervasive 
in the quantum hydrodynamics of incompressible and com- 
pressible quantum Hall systems, as a direct expression of their 
quantum hydrodynamics Similarly, local gauge symme- 
tries appear naturally in the context of strongly correlated sys- 
tems such as the t — J model, quantum dimer models, and 
other systems |7|]. Of particular interest for the problems dis- 
cussed in this paper are the sliding phases of arrays of Lut- 
tinger liquids @|, quantum Hall smectic (stripe) phases 0], 
DNA intercalates in lipid bilayers lHoll . as well as in some ring 
exchange models of frustrated antiferromagnets Fill . The dis- 
crete sliding symmetries we discuss here are a discrete, Z2, 
version of the continuous sliding symmetries of the systems 
mentioned above. The existence of sliding symmetries has 
profound effects on their quantum phase transitions ? whose 
behavior only begun to be understood quite recently IIFZll and 
still remains largely unexplored. 

Amon gst o thers, discrete sliding symmetries are present in 
spin OullH, orbital_|SEl^§[i§|23] 5 an d supercon- 
ducting array systems 12ll UM . We further demonstrate that 
the planar orbital compass model and the Xu-Moore 

model of two dimensional p + ip superconducting 

arrays are, in fact, one and the same system, related by a sim- 
ple duality transformation. Viewed in that light, the discrete 
sliding symmetries which the Hamiltonians describing super- 
conducting arrays display are natural. By applying our new 
dualities, we find self-dualities for the three dimensional or- 
bital compass model and several other systems. These du- 
alities do not rely on operator representations ll23ll . l24ll nor 
on standard combinatorial loop/bond counting arguments or 
summation formulas 12511 . Rather, the new dualities that we 
report here appear as simple geometrical reflections between 
various spin and spatial axis. The dualities investigated in 
this paper map such trivial geometrical reflection self-dualities 
in one model onto far less trivial weak-strong coupling self- 
dualities in other systems. In a formal setting, our dualities 
correspond to different space-time cuts of a single classical 
action. Choosing a certain time and spin quantization axis, 
we find one spatial system whilst choosing the time axis to 
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lie along another direction in space-time leads to a seemingly 
very different (yet dual) spatial model. Our new, purely geo- 
metric, dualities further extend and complement, from a rather 
general perspective, the dualities generally derived via tech- 
niques such as those in, e.g., |2(J. 

The plan to the paper is as follows. In Section[n] we intro- 
duce the planar orbital compass model in both its isotropic and 
anisotropic incarnations. We identify the many gauge like and 
single reflection symmetries of this model (the latter reflec- 
tion symmetry will, as we will later find out, play the role of 
self-duality). In the aftermath, we construct order parameters 
invariant under these symmetries. 

In Section[[!3] we discuss another two dimensional XY sys- 
tem that possesses one dimensional gauge like symmetries. 
This system has been argued to embody the quintessential 
physics of a square lattice array of p + ip superconducting 
grains. As we will show, this model is identical to the planar 
orbital compass model. 

In Section IIVI we discuss the three dimensional orbital 
compass system. This model has been considered to embody 
the prototypical features of orbital system Hamiltonian and 
might be directly relevant to the so-called "£29" systems (such 
as the vanadates and manganates) in particular. We identify 
gauge like symmetries in this system. As in the planar case, 
we find nematic orders invariant the gauge like symmetries. 
This in turn suggests that orbital systems might possess ob- 
servable nematic orders. 

In Section|V] we employ simple geometrical reflections to 
derive dualities for extended systems (now resi ding in three 
dimensions). The dualities of Xu and Moore fell 12211 (derived 
by Kramers-Wannier loop counting) form a subset of the de- 
rived dualities. The central actor in our scheme is a geomet- 
rical inversion operator which allows us to set the imaginary 
time axis along different external space-time directions with 
similar ideas for choosing the internal spin quantization axis. 
These operations generate, in turn, many different dual mod- 
els. With this geometrical understanding of the observed du- 
ality in hand, we return to the self-dual point of I2H I22I1 and 
make comparisons to other systems. 

In Section lvTI we use an operator representation of the dual- 
ity transformation to rederive the dualities for the orbital com- 
pass model which we obtained via geometrical reflections in 
the previous sections. Section IVIll is devoted to the conclu- 
sions. In Appendix IaI we discuss the self-duality of "around 
the cube" models in transverse fields. 



II. QUANTUM PLANAR ORBITAL-COMPASS MODELS : 
SYMMETRIES. 

We start with the planar compass-model. The compass 
models often serve as the simplest caricatures for the physics 
of 3d orbital systems wherein Jahn-Teller interactions as well 
as magnetic exchange processes are dictated by the orienta- 
tion of the orbitals at the various lattice sites. In the orbital 
compass models, the spin variables code for the orbital states. 
As orbitals extend in real space, all orbital dependent inter- 
actions are highly anisotropic- these interactions link the ex- 



ternal lattice directions with the internal "spin"_fie- orbital) 
orientations. We refer the interested reader to 12711 where the 
physics of orbital systems and the orbital only models that we 
investigate is explored in depth. 

The planar compass model is defined on the square lattice 
where at each site r there is a S= 1/2 operator denoted by S r = 
7rcr r . The isotropic planar orbital model Hamiltonian 

H iso = -J^2K<T$ +ex +a*a* + J, (2.1) 

r 

where the nature of the interaction allows us to set J > 1 28 ] . 

Unlike the more conventional nearest neighbor spin Hamil- 
tonians which posses a continuous global rotational symme- 
try, the compass model Hamiltonian is not invariant under ar- 
bitrary global rotations of all spins. [Physically, the lack of 
this symmetry is the direct consequence of the coupling be- 
tween the internal polarization directions (orbital states) and 
the external lattice directions (as much unlike spins, the or- 
bitals extend in real space) 12711 .1 Instead, this model pos- 
sesses many new non-trivial symmetries corresponding to 
specific quantized angles of rotation of all spins on given 
rows/columns and a single additional rather trivial reflection 
symmetry (which upon mapping will enable us to find a non- 
trivial weak to strong coupling self-duality in another model). 
As a consequence of these symmetries, this model harbors an 
infinite degeneracy of all states and of its ground states in par- 
ticular. Let us consider the system with open boundary condi- 
tions on an L x L square lattice and let us define an operator on 
an arbitrary horizontal line (of ordinate z) O z = Ylx=-L a x 
and an operator on an arbitrary vertical line (of horizontal 

intercept x) O x = YI^=-l <j z- ^ * s readily verified that 
for all sites r whose z component is r z = z, the product 
6j 1 a56 z = -<j5 while OjVSO, = 0%. Similarly, 6 X 
inverts the z component of all spins on a vertical line, while 
leaving a x untouched. In the case of symmetric exchange 
constants for bonds along the x and z axis, as in the com- 
pass model under consideration here, (where both exchange 
constants are equal to J), we further have a single additional 
Z2 reflection symmetry (<j x — > <r z , a z — ► cr x )- a rotation by it 
about the symmetric line (the 45 degree line in the xz plane), 

i.e. 6 Reflection = Ilf ex P [* ^ (°f + °*0] ■ For each of these 
operations, O^HOa = H. This symmetry, 6 Reflection is 
a manifestation of self-duality present in the model- we will 
explain the origin of this comment later. Putting all of the 
pieces together, as a consequence of these symmetries, each 
state is, at least, C(2 L ) degenerate. Formally, these symme- 
tries constitute a gauge like symmetry which is intermediate 
between a local gauge symmetry (whose volume scales as the 
system area) and a global gauge symmetry (whose logarithmic 
volume is point-like). These intermediate gauge symmetries 
suggest that non-trivialities may occur. As it turns out, such 
large discrete symmetries do not prohibit ordering in classical 
variants of this model albeit complicating matters significantly 
fl7l [I9I.I20I1 . This ordering tendency may be expected to be- 
come fortified by quantum fluctuations ("quantum order out 
of disorder") 12911 . 

As an aside, we note that the global nematic symmetries 
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(the global rotation of all spins by nw/2 with n = 0, 1, 2, 3) 
are not independent symmetries on top of the gauge like sym- 
metries discussed above. Rather, there is only one (<r x — > a z , 
a, — » a x )- a rotation by it about the symmetric line (the 45 
degree line in the xz plane)) additional symmetry supplant- 
ing the gauge like symmetries. To see this, first note that the 
global inversion operation, a — > —a, is a composite of the 
row/column inversion symmetries: An inversion of a x on all 
rows followed by an inversion of a z on all columns leads to 
the global inversion operation. Next, note that by fusing the 
global inversion symmetry with the global Z 2 reflection sym- 
metry, we may produce the four global nematic symmetry op- 
erations (rotations by mr/2). Thus, unlike what is suggested 
by l20ll . the global nematic symmetries do not supplant the 
gauge-like symmetries and no less important, the quantum 
system possesses the above reported gauge like symmetries 
embodied by the operators X)Z . 

The classical (large S) ground state sector of the orbital 
compass model further possesses an additional continuous 
(C/(l)) symmetry not captured by the discrete (Z 2 ) 2L+1 sym- 
metries (2 L of these associated with horizontal, 2 L associated 
with vertical discrete spin flip symmetries, and one Z 2 sym- 
metry being the (a x — > a z , a z — > a x ) reflection symmetry) 
detailed above. This continuous symmetry is made evident by 
noting that any constant spin-field, a r = a, is a ground state. 

First, we note that 2~2 a=x »[ <7 r ] 2 is constant. Thus, up to an 
irrelevant constant, the general Hamiltonian of Eq. ( 12. Q is 

VAX 

which is obviously minimized when the spin field is constant. 
We emphasize that the continuous symmetries which under- 
score these ground states are just symmetries of the states and 
not of the Hamiltonian itself. In common parlance, these are 
emergent symmetries specific only to the ground state sector. 
Therefore, at least in the orbital-only models, we are not in a 
setting where a Mermin- Wagner argument can be applied. 

With an eye towards things to come, let us now introduce 
and examine the anisotropic planar compass model, 

r r 

It is readily verified that this more general Hamiltonian har- 
bors all of the one dimensional gauge like symmetries encap- 
sulated by O x . z . The only symmetry which does not persist 
for arbitrary J x , J z is the reflection symmetry. [Insofar as 
its underlying physics is concerned, this anisotropic compass 
model emulates Jahn-Teller distortions on a strained lattice 

The two terms in the anisotropic compass model of 
Eq. J2.3l >. trivially compete. The first term favors ordering of 
the spins parallel to the x axis while the second favors an 
ordering of the spins parallel to the z-axis. Order becomes 
more inhibited when the competition between the two terms 
becomes the strongest (J x = ± J z ) as it indeed occurs within 
the compass model of Eg. ( 12. 11 1. We note that the gauge like 



symmetries (encapsulated by the column/row O x>z genera- 
tors) preserve the Hamiltonian also for arbitrary \ J X \ ^ \J Z \. 
A natural (smectic-like) order parameter in the orbital com- 
pass model monitors the tendency of the spins to order along 
their preferred directions 

(Just as in smectic liquid crystals, having all spins point in 
the e Q direction or in the (— e a ) direction is one and the same 
insofar as the above order parameter is concerned). This ne- 
matic like order parameter is invariant under all gauge like 
symmetries. 

Similar to the xy symmetric order parameter above, for the 
anisotropic planar orbital compass model (say \J X \ > | J z \), 
we may consider the Ising like nematic order parameter, 

m x = (a^ +g J, (2.5) 

with a similar definition for the system with \J Z \ > \J X \. 
These order parameters are invariant under the gauge-like 
symmetries of the system. 

The classical, large S, rendition of this model, has simi- 
lar nematic like order parameters invariant under all gauge- 
like symmetries fl7t fl^. Eoll . Here, and in fact for all spins 
S > 1/2, the order parameter can be local (not a bond or- 
der parameter involving two spins). All quantities Q a p = 
gag/3 _ LS a p, with a, (3 internal spin indices, and with d 
the dimension (d = 2 in the planar higher spin extensions 
of the orbital compass model) are invariant under all gauge 
like symmetries. The order parameter (Qu) is anticipated for 
\Jx\ > \Jz\ (and (Q22) for \ J Z \ > \ J X \). Similar quantities 
may be introduced for higher dimensional (d > 2) generaliza- 
tions of the planar compass model. 



TIL p + ip SUPERCONDUCTING ARRAYS 

A Hamiltonian describing a square lattice of p + ip su- 
perconducting grains (e.g. Sr 2 Ru04) was recently suggested 
EH, 

H = -K^2<T z cr z a z a z -h^a*. (3.1) 

□ r 

Here, the four spin product is the product of all spins residing 
at the four vertices of a given plaquette □ (not on its bonds 
as for gauge fields!). As noted by Xu and Moore, l2lll . the 
quantity 

Op = n CT - (3-2) 

r 

with the string product (along "P") extending over all spins in 
a given row (r z = z) or a given column (r x = x), is con- 
served. The discrete (gauge-like) sliding symmetry of this 
model is similar to that of the planar orbital compass model 
and we will indeed show that these two models are actually 
dual to each other. 
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The central derivation in I2H I22I1 was a self-duality of the 
Hamiltonian in Ea. d3.U via a tour de force Wannier Kramers 
loop counting arguments. The form of this self-duality is 
somewhat similar (yet still very different) to the beautiful self- 
dualities of l3lll . Similar dualties were discussed in the ring 
exchange systems of full . In these models not only a relation 
amongst strong and weak coupling is given by the self-duality 
but the self-duality further intertwines the various terms (e.g. 
large h is related to small K in the self-duality of Eq. (I3.1> and 
vice versa as found by Xu and Moore). 

We will shortly establish that the rather complicated look- 
ing weak coupling to strong coupling self-duality of Ea. (l3.1> 
derived by Il2ll 12211 immediately follows from a very simple 
purely geometric (Z2 reflection) self-duality of the planar or- 
bital compass model. This self-duality may also be related 
(albeit in a less general fashion) to the trivial geometrical self- 
duality of the planar orbital compass model via the operator 
representations of Section dVll . In the aftermath, the plaquette 
coefficient K in Eq. d3.H may be related to the exchange am- 
plitude J x of Eg. (12. 3> whereas the transverse magnetic field h 
of Eq.( l3.1> becomes trivially related to J, of Eq. d2.3> . 

IV. SYMMETRIES OF THE THREE DIMENSIONAL 
ORBITAL COMPASS MODEL 

The canonical prototype of all orbital-spin 0] and orbital- 
orbital interactions is the orbital compass model l30fl . The 
model is defined on the cubic lattice where at each site r there 
is an S=l/2 operator denoted by Sp = The orbital model 
Hamiltonian 

H = j]T(a5^ + olo%^ + <J^ + J. (4.1) 

r 

Let us define an operator on an arbitrary xy plane P (of in- 
tercept z) Op- z = IlreP a ? w i m simile definitions for Op, x 
and Op vy . These operators may be recast as rotations by it 
about an axis orthogonal to the plane. For all sites f in the 
xy plane P whose z component is r z = z, up to a multiplica- 
tive phase factor, the operator Op ;z = exp[i(Tr/2)crp/H] with 
a p = J2reP °F- The P roducts 0~ P ] z a^ ' x ' v d P . z = -ap y 
while Op^jjpOp-z = <7p. Similarly, Op- x inverts the y and 
z component of all spins on the yz plane of intercept x while 
leaving a x untouched. These "string" operators spanning the 
entire plane commute with the Hamiltonian, [H, Op- a ] = 0. 
The classical orbital compass model has an exact [Z2] 3L 
symmetry (along each chain parallel to the cubic a (x, y, or 
z) axis, we may reflect the a spin component, S a — > — S a , 
while keeping all other spin components unchanged, Sp^ a 
S/3^ a )- The quantum orbital compass model has a lower exact 
[Zy 3L gauge like symmetry (forming a subset of the larger 
py 31 symmetry present for classical spins). As alluded 
to above, the gauge like [Zy 3i symmetries of this quantum 
S = 1/2 case (as well as all representations), become evident 
once we rotate, with no change ensuing in the Hamiltonian, 
all spins in a plane orthogonal to the cubic lattice direction a 
by 7r about the internal S a quantization axis. 



As before, let us now introduce and examine the anisotropic 
orbital compass model, 

H = ~Jx^2 a T a r+&* ~ J V ^2 a r CT r+e y ~ J z X! ^r+e* • 
r r r 

(4.2) 

(The isotropic orbital compass model corresponds to J xy z = 
— J). The anisotropic orbital compass model possesses all of 
the gauge like symmetries of the isotropic orbital compass 
model (planar rotations in the quantum model and more nu- 
merous single line inversions in the classical case). Further, 
if at least any two of the three exchange constants { J a } are 
identical the system possess a reflection symmetry. 

As in the planar orbital compass model, nematic like or- 
der parameters may be constructed for both the isotropic and 
anisotropic systems. Thus, we naturally predict the existence 
of observable nematic orbital orders in t2 g systems. 

In what follows, we will also investigate a related system 
governed by the Hamiltonian 

H = -J x ^2 Cr ? Cr ?+e x - J V ^2 a i a r+e y ~ -4 / ] CT r ^r+e, ■ 
r r r 

(4.3) 

Note the similarity between the XY model of Ea. d4.3t and 
the orbital compass model of Ea. d4.2t . In the limit J z = 0, 
Ea. J4.3b trivially degenerates into the strained planar orbital 
compass model of Eq. d2.3> . 

We will construct new "plaquette models" (in which the 
spins reside on the lattice sites not on bonds) dual to Ea. d4.3l > 
which possess a self-duality and gauge like symmetry, natu- 
rally extending the results of I21II22I1 . 



V. NEW DUALITIES AND SELF-DUALITIES FOUND BY 
PLANAR REFLECTIONS 

We now transform the zero temperature Quantum problem 
of EQ. dXTt onto a classical problemin (d+ 1) dimensions. To 
this end, we work in a basis quantized along cr 2 (= ±1). We 
now consider the basis spanned by two spins (cr 2 , a zl ) at the 
same spatial site r yet at two consecutive imaginary times t 
and (r + At). The transfer matrices corresponding to ae ha 
(stemming, in the imaginary time formalism from a propaga- 
tor e~ HAT such as e ' lffXAr with h = hAr) and e 1 ' 7 *' 7 " (or, 
with space time coordinates explicitly instated, e J<Jl - T<Jl - T +* T ) 
are the same provided that tanh/i = e~ 2J (or equivalently 
sinh27isinh2J = 1) and a = (2sinhJ) 1/2 . Similarly, the 
non-vanishing eigenvalues of the transfer matrices 

CXp[/i,<X+l,r<T+AT CT m,r] (5-1) 

and 

exp{J x al T (7f +ltT ] (5.2) 
are equivalent once sinh 2K X sinh 2 J x = 1, 
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FIG. 1 : The classical Euclidean action corresponding to the Hamil- 
tonian of Eq, <3.1> at zero temperature in a basis quantized along the 
<j z direction. The transverse field leads to bonds parallel to the imag- 
inary time axis while the four plaquette interactions become repli- 
cated along the imaginary time axis. Taking an equal time slice 
of this system we find the four spin term of Eq. J3.1> and the on- 
site magnetic field term. If we interchange r with z, we find the 
anisotropic planar orbital compass model of Eo. <2.3> in the basis 
quantized along the <j x direction. 



In the standard imaginary time mapping of quantum sys- 
tems to classical actions, we identify J a = J a Ar with the 
aforementioned At the lattice spacing along the imaginary 
time direction. 

The generalized classical Euclidean action corresponding 
to Eq. d2.U is 



C JS \ ^ z z z z 

° ~ lx x CT r,T Cr r,T+AT r+e 3c ,T°r+e X: T+AT 

□ e(a:T) plane 

r 

(5.3) 



A schematic of this action in Euclidean space-time is shown 
in Fig.Q. If we relabel the axes and replace the spatial in- 
dex x with the temporal index r, we will immediately find 
the classical action corresponding to the the Hamiltonian of 
Eq. d3.ll > depicting p + ip superconducting grains in a square 
grid. This trivially suggests that the anisotropic planar orbital 
compass system (Eq. d2.3» and the Xu-Moore Hamiltonian 
(Eq. d3.ll l) are dual to each other. In Section dVU . we sketch 
a detailed derivation of this duality by the operator dualities 
of l23ll . l24ll . This classical action follows from the equiva- 
lence of the transfer matrices corresponding to Eqs. d5.1l5^ 
or, alternatively, from the equivalence of Eq. d2.3l l) to Eq. d3.1> 
[which will be proved in detail by operator representations 
in Section fVll il and the relation between the transfer matrices 
corresponding to e ha and e ,1a a ~ . 

We find that the classical action corresponding to the model 



ofEq.O is 

S=[-tan\T 1 {e- 2J ^ T ) ^ aaaa 
□eir plane 

— AtJ z aa 
z direction 
-tanh-^e" 2 ^^) ^ aaaa}. 

Ueyr plane 



(5.4) 



Here and elsewhere, a = ±1 are c-numbers and we omit the 
(z) polarization superscripts. 

We now extend the duality of self-duality of Eq. O.H to the 
three dimensional arena. First note that by interchanging the 
imaginary time coordinate r with the spatial z coordinate, we 
find that 



(5.5) 



H = -{K xz Y, v z cr z a z a z 
-K yz J2 v z v z <? z v z + hJ2<) 

Deyz r 



is dual to the system given by the Hamiltonian of Eq. d4.3> . Let 
us now derive self-dualities of this extended three dimensional 
system (en passant, effortlessly proving the central result of 

El). 

Expressing the action corresponding to the Hamiltonian of 
Eq.( l4^ in a spin eigen-basis of a x and inverting the spatial z 
and x coordinates of any site r (a Z2 operation), we obtain 



Sduai = [~ tanh (e 



□giT plane 

—AtJ x aa 
z direction 
-tanh- 1 ^- 2 ' 7 ^) Y aaaa\. 

□eyr plane 

Looking at Eqs.( l5.4l5.6> , we see that if 

S = [A aaaa 
□ en plane 

+B Y aa 
z direction 

+C aaaa] 



(5.6) 



(5.7) 



then 



□ej/T plane 

Sdual =[A Y aCT(Ja 

Oexr plane 

+B Y aa 
z direction 
+C aaaa}. 
Oeyr plane 



(5.8) 
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Here, A = - tanh -1 ^- 2 - 7 *), B = -J z , and C = 
-tanh" 1 (e" 27 «). Similarly, A = -tarm -1 ^ -27 *), B = 
— J x , and C = C. Eliminating J x . z , we find that 

sinh2isinh2B = 1, 
sinh274sinh2B = 1, 

C = C. (5.9) 
Taken together, these relations imply that 

sinh2Asinh2B = 1 (5.10) 

is a self-dual line for any value of C. 

This extends the dualities of I21H22I1 in a very natural fash- 
ion to higher dimensions. Moreover, note in this formalism 
the dualities just "fall into our lap"- no involved calculations 
nor loop counting were necessary. The duality is a trivial ge- 
ometrical reflection. 



VI. OPERATOR DUALITY TRANSFORMATIONS OF THE 
QUANTUM PLANAR ORBITAL-COMPASS MODELS ONTO 
QUANTUM ISING PLAQUETTE MODELS 

A duality between the planar orbital compass model 
(~Ea. J2.3l l) and the Xu-Moore model of Eq.l D.H is suggested 
by the cubic point group operation interchanging x with r in 
Fig.Q). We now prove this duality at all temperatures. To 
this end, we invoke a simple operator duality transformation 
followed by a summation over the horizontal bonds (which 
amounts to a standard gauge fix) in the model that results. 
The upshot of the up and coming discussion is that the quan- 
tum planar compass model of Eqs.l l2.ll i can be mapped onto 
the Hamiltonian of Eq.l l3.ll l precisely at its zero temperature 
self-dual point. 

The salient feature of the Pauli matrices a x and a z is that 
they anti-commute at a common site while commuting every- 
where else. It is readily verified l23Tl . l24ll that these relations 
are preserved by the canonical duality relations on the dual 
lattice 

— Z ^X X X X f£l 1\ 

<7 T = T T T T (0.1) 

with the plaquette product of t x on the right hand side corre- 
sponding to the four spins surrounding the dual lattice site r* 
(the center of the plaquette as shown in Fig.© below), and 



n 

x<x* 



x,x-\-e z 1 



(6.2) 



the product of t z placed along vertical bonds (linking x and 
x+e z ) along a horizontal line- see Fig.(|3}. The series of trans- 
formations below leading to Eq.( I6.4> may be vividly followed 
in Figures ^E) . 

Inserting Eqs. j6. 116.21 1 into Eq. j2.1> . we obtain 



Hi, 



-/Vfr 2 



7 , r r*,r*+e z T r*+2e z ,r*+2e z T r*+2e z ,r*+2e z 



he z — e x ,r* — e x ^~r* — e x ,r* 

















* 



























FIG. 2: The lattice and a dual lattice site (marked by an asterisk "*" 
at a plaquette center). Here we illustrate the representation of o z as 
on a given dual plaquette site as the product of r x operators placed 
on all 4 bonds composing the plaquette of the original lattice. 



The first term corresponds to an external transverse magnetic 
field of strength J along the z-axis acting on all vertical bonds 
while the second term encapsulates the product of six bonds 
forming the outer shell of two plaquettes pasted together along 
the z-axis. The bond common to the two plaquettes evapo- 
rated due to the relation r 2 = 1. The net result of Eq. d6.3l l is 
shown in Fig.Q. 

Next, we choose the longitudinal gauge wherein all hor- 
izontal bonds have t z = 1. This can be achieved via ex- 
plicit gauge transformations or by simply noting that in the 
representation with horizontal bonds with t z = 1, the dual- 
ity relations of Eqs. J6. II I6.2l i become identical to the dual- 
ity relations in one dimensional spin chains (performed inde- 
pendently for each horizontal row) which trivially satisfy the 
commutation of spin variables of different sites, the anticom- 
mutation of the x and z components of the spin on the same 
site and the square of each spin operator. In this longitudi- 



nal gauge where if. 



r*+e x 



1, the Hamiltonian now involves 



only vertically oriented bonds (parallel to the z-axis). Defin- 



ing spins 



H. 



r* +e z r* +e z — e x r" e x 



-hJ24*, (6-4) 



(6.3) 



with the new parameters h = K being equal to the former J 
of Eq. J6.3> . Thus the isotropic planar compass orbital model 
lies precisely on the zero temperature self-dual line ft, = fc of 
Eq.( l6.4l . This result is shown in Fig.(|5}. 



VII. CONCLUSIONS 

In conclusion, we investigated several systems displaying 
discrete (gauge-like) sliding symmetries and illustrated that 
two such systems are dual to each other. The enhanced dis- 
crete sliding symmetries in these systems go hand in hand 
with a dimensional reduction that occurs in several limiting 
cases of these systems (e.g., the system of Eq. J3.lt in the limit 
of ft = is none other than 1+1 dimensional version of the 
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* 



























FIG. 3: A graphical representation of <j x the string product of r z on 
all vertical bonds from the boundary up to the dual lattice site. 

















* 



























FIG. 4: The product of Jx<J x a x +6x becomes J x t z on a single verti- 
cal bond on the right-hand side of the dual plaquette site correspond- 
ing to r. The product Jz^r^r+e^ becomes in the dual representation 
the product of all r x operators forming the outer shell of a verti- 
cal domino multiplied by J z . Putting all of the pieces together, the 
Hamiltonian becomes the sum of J z multiplying a domino shell of 
t x on bonds augmented by J x multiplying a single vertical bond on 
which t z is placed. 





i 


> . 








t 





























FIG 5: Choosing a gauge in which r x — 1 on all horizontal bonds, 
identifying the centers of the vertical bonds as sites, we find that the 
Hamiltonian corresponds to the product of four tj x operators on the 
vertices of a plaquette [K s x * s x * +g z s?* +e z -e x s r* -S x ] augmented 
by an external transverse field giving rise to hs z , . Here, h = J x 
and K — J y . Thus, the planar orbital compass model is dual to the 
superconducting array system of I2TH22I1 . 



that this is an unconventional quantum phase transition whose 
universal behavior is worth understanding. 
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one dimensional Ising model.) We found nematic order pa- 
rameters invariant under these symmetries. This suggests the 
specter of detectable orbital nematic orders in t2 g orbital sys- 
tems. 

The superconducting array model of I2II1 . 122I1 is dual to the 
planar orbital compass model and as such has a finite temper- 
ature transition for a large S incarnation at its self-dual point. 

We find that dual models may be derived by flipping the 
spatial and imaginary time axis (and/or quantization axis). In 
an upcoming work, we will elaborate on this novel approach 
to dual models as different cuts of a higher dimensional theory 

The nature of the quantum phase transitions in these sys- 
tems remains an open problem. A straightforward examina- 
tion of the Xu-Moore model shows that the finite temperature 
transition from a high temperature disordered phase to a low 
temperature phase in which the product of Ising spins on pairs 
of sites belonging to different sub-lattices orders. This classi- 
cal transition is continuous, and in the universality class of the 
2D classical Ising model. Although the nature of the T = 
transition is not established at present time, there are sugges- 
tions that it may actually be a continuous quantum phase tran- 
sition. The presence of the discrete sliding symmetry suggests 



APPENDIX A: SELF-DUALITY OF "AROUND THE CUBE" 
MODELS IN THE PRESENCE OF TRANSVERSE FIELDS 

In this appendix, we explicitly generalize the self-duality 
that we obtained earlier for plaquette models with a transverse 
field to cubic models with eight spin interactions au gme nted 
by a transverse field. Such a duality was alluded to in I2ll 12211 . 

In the below, we derive this duality by going back and forth 
from various quantum systems to corresponding (d + 1) di- 
mensional classical actions when different spin quantization 
and spatial lattice directions are chosen. 

To prove the self-duality of such cubic systems, first con- 
sider the Hamiltonian 

H = -K o v o v o y o y 
□ emplane 

—Jz J2 (J Z (J Z . (Al) 

bonds along z axis 

If we write down the classical action in a spin basis quantized 
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along the o z axis, we find 

S\ = — A'i oooooooo 
cubes in xyr 

-Jzi Yl ua - (A2) 

bonds along z axis 

Here sinh 2K i sinh 2K = 1 with K = KAt, and J z \ — 
J z At. Alternatively, if we write down the classical action 
corresponding to Eq. JAH in a spin basis polarized along o y , 
we find 

S 2 = —K'2 oooo 
□exyplane 

-Jz2 ^ oooo. (A3) 
□ezrplane 

Here, K2 = KAt and sinh2J Z 2 sinh2J z = 1 with J z = 
J z At. Thus, we find that the classical actions Si and 5*2 are 
dual to each other. The classical action Si also corresponds to 
the Hamiltonian 

TT Ts- \ ^ „Z Z Z Z Z Z Z_Z 

"cube — / OOOOOOOO 

cubes 111 xyz 



when written in a spin basis quantized along o z . Thus, H cu b e 
may be represented by the classical action S 2 - Putting all of 
the pieces together we find that 

sinh 2if* At sinh 2if 2 = 1 
smh2h z *Ar sinh2J z2 = 1. (A5) 

Interchanging, in the action 52, K — * J, x — ► z, y — * 
t, we obtain a new action (S3) whose partition function is 
identically the same. By the same steps outlined above, the 
classical action S3 corresponds (via duality transformations) 
to the Hamiltonian 

TJ IS \ ^ „z z z z z z z z 

H-cube = ~ -"■* / OOOOOOOO 

cubes in xyz 

r 

This establishes the duality between H cube and H cu b e , 
fu = h„ 

h z , = K«. (A7) 

Fusing these relations together, we find that ft* = K„ consti- 
tutes a self-dual line of H cu b e . 
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